NONLINEAR SCHRODINGER EQUATION AND FREQUENCY SATURATION 
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Abstract. We propose an approach that permits to avoid instabihty phenomena for the 
nonlinear Schrodinger equations. We show that by approximating the solution in a suitable 
way, relying on a frequency cut-off, global well-posedness is obtained in any Sobolev space 
with nonnegative regularity. The error between the exact solution and its approximation 
can be measured according to the regularity of the exact solution, with different accuracy 
according to the cases considered. 



1. Introduction 

We consider the nonlinear Schrodinger equation 

(1.1) z9tu + Alt = e|up'^u, (i, x) e / X R''; U|t=o = "o, 

for some time interval I 3 0, with e — 1 (defocusing case) or e = —1 (focusing case). The 
aim of this paper is to propose an approach to overcome the lack of local well-posedness 
in Sobolev spaces with nonnegative regularity. 

Recall two important invariances associated to ( II. Il l: 

• Scaling: if u solves (II. lb . then for A > 0, so does ux{t, x) :— \^/"'u (A^i, Ax). 
This scaling leaves the iJ|=-norm invariant, with = d/2 — 1/cr. 

• GaHlean: if u solves ( II. Il l, then for v e R"^, so does e™ '^~*l''l^*/^u(i, x — vt). 
This transform leaves the L^-norm invariant. 

These two arguments suggest that the critical Sobolev regularity to solve ( II. Il l is max(sc, 0). 
Indeed, if Sc ^ 0, local well-posedness from H''{IV^) to H^{Il'^) for s ^ Sc has been es- 
tablished in |13|, and if Sc < 0, local well-posedness from H''(R'^) to H''(R'^) for s ^ 
has been established in ||27l . 

If Sc > 0, pathological phenomena have been exhibited for initial data in H''{'El'^) with 
< s < S(,'- Gilles Lebeau has proved a "norm inflation" phenomenon for the wave 
equation dfu - Au + = 0, x e R^, p g 2N + 1, p ^ 7 (IHI; see also [24]). The 
analogous result for (II. Il l has been established in lfT4l and H. 

Theorem 1.1 (From pi, 041). Let cr ^ 1. Assume that Sc = d/2 - 1/a > 0, and let 
< s < Sc- There exists a family (mo)o</i<i in SiW^) with 

llwollff^K") Oas /i 0, 
a solution m'' to ( II. lb and < ^'^ — ?> 0, such that: 

-> +00 as h ^ 0. 
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The argument of the proof consists in considering concentrated initial data, 




and showing that for very short time, the Laplacian can be neglected in ( 11.11 ). The above 
result then stems from its (easy) counterpart in the ODE case, by choosing a suitable a > 0. 
In the spirit of |23|, the above result has been strengthened to a "loss of regularity" in 
lUllHll^l; the assumptions and conclusion are similar to that in Theorem ll.il the only 
difference is that u^{t^, ■) is measured in iJ'^(R'^) for any k > s/(l + a{sc — s)), thus 
allowed to be smaller than s. In all the cases mentioned here, the lack of uniform continuity 
of the nonlinear flow map near the origin is due to the appearance of higher and higher 
frequencies on a very short time scale. If Sc < 0, similar pathological phenomena have 
been established in H'^ (R'') with s < 0, where on the contrary, low frequencies are ignited; 
see e.g. |l2l[T0l[l4]|20l. In the rest of this paper, we focus on nonnegative regularity, s ^ 0. 

The goal of this paper is twofold. First, we want to investigate a method to remove the 
pathology mentioned above, causing a lack of well-posedness for dl.U . in a deterministic 
way, as opposed to the probabiUstic approach initiated in ||6l|2l for the wave equation. The 
other motivation is related to numerical simulations for ( II. Il l, where high frequencies may 
be a source of important errors; see for instance [18,1 , a reference which wiU be discussed 
further into details in Sections[3]and|4l 

We show that with a suitable cut-off on the high frequencies of the nonlinearity, the ob- 
structions to local well-posedness vanish, and the problem becomes globally well-posed: 
the nonlinear evolution of any initial datum in L^(R'^) can be controlled a priori, an in- 
formation which may be useful for numerics, since we do not have to decide if the initial 
datum belongs to a full measure set or not. This strategy is validated inasmuch as this 
procedure yields a good approximation of the solution to (II. lb as the cut-off tends to the 
identity. Note that this approach can be viewed as a deterministic counterpart of the one 
presented in [5 | (see also |3|). There, for the one-dimensional -supercritical defocusing 
nonlinear Schrodinger equation, the authors construct a Gibbs measure such that, among 
other features, the pathological phenomenon described in Theorem 11.11 occurs for a set 
of initial data whose measure is zero: on the support of the Gibbs measure, the Cauchy 
problem is globally well-posed, and a scattering theory is available. Both points of view 
aim at showing that norm inflation in the sense of Theorem ll.ll is a rare phenomenon: in 
||5], the authors give a rigorous meaning to this statement in an abstract way, while we are 
rather interested in a recipe to avoid instabilities for sure, by a suitable approximation of 
the equation, which can be used typically for numerical simulations. 

Our choice of cutting off the high frequencies instead of, for instance, the values of 
the function itself is indeed motivated by numerics, where it is standard to filter out high 
frequencies (sometimes without even saying so). In an appendix, we discuss another ap- 
proach, consisting in saturating the values of the nonlinearity. One could of course combine 
both approaches, frequency and physical saturations. 

Notations. We define the Fourier transform by the formula 



We write a < 6 if there exists C such that a ^ Cb. In the presence of a small parameter h, 
the notation indicates that C is independent of ft. £ (0, 1]. 
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2. From instability to global well-posedness 

Let X '■ R-'' ~^ [0, 1] be a smooth function, equal to one on the unit ball, and even: 
x{~x) = x{^) for ^11 X G R*^. It may be compactly supported, in the Schwartz class 
5(R'^), or with a slower decay at infinity. For simplicity, we will not discuss sharp as- 
sumptions on X- We define the frequency "cut-off" 11 as the Fourier multiplier 

As pointed out in the introduction, in the examples constructed to prove the lack of local 
well-posedness, the mechanism of high frequencies amplification occurs at the level of the 
ordinary differential equation. We discuss some strategies to saturate high frequencies at 
the ODE level first, with e — 1 for simplicity. 

2.1. Candidates at the ODE level. The first possibility to prevent the appearance of high 
frequencies by nonlinear self-interaction consists in saturating the whole nonlinearity: 

(2.1) idtv = Ti{\v\^''v) . 

This can be viewed as an extremely simplified version of the /-method (see e.g. fTSl). 
Another choice consists in saturating the high frequencies in the "nonlinear multiplicative 
potential" only, that is for cr £ N, we propose two possibilities, 

(2.2) idtv = n V, 

(2.3) idtV^{\l{\v\^)Y V. 

In the cubic case (7=1, the last two approaches obviously coincide. These two approaches 
have two advantages over ( 12.1b : 

• They preserve the gauge invariance. If v solves the equation, then so does ve^^ for 
any constant e R. 

• They preserve the conservation of mass. 

To see the second point, rewrite n(/) = K * f, with K{x) = (27r)^''/^x(^2;). Since x 
is even and real-valued, so is K, and therefore 9t|wp = in ( 12.21 ) and ( I2.3l l. This identity 
leads to the conservation of the L^-norm at the PDF level. 

Before passing to the PDF case, we conclude this section by showing that even at the 
ODF level, cutting off high frequencies in the initial data does not suffice to prevent the 
appearance of higher frequencies in the solution for positive time. For a G iS(R'') and 
s > 0, consider the solution to 

Then I'lJ^g is bounded in H^(R!^), uniformly in h G (0,1], and if a is compactly supported 

(in i?(0, R)), then u''|t=o is compactly supported (in i?(0, R/h)). Since 9f — 0, we 
have the explicit formula 

v''{t,x) = h'-^'^a {£j cxp (^-zi;i2.(.-<i/2) 
We check that for t > 0, as /i — s> 0, the homogeneous Sobolev norms behave like 

at least for fc G N. The above quantity is unbounded as /i — ^ if 



l + 2CT(s-d/2)' 
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Therefore, if s < d/2, v'^{t, •) is unbounded in iJ*(R'^) for t > 0, as h ^ 0: cutting 
off the high frequencies in the initial data does not suffice to control the frequency support 
of the solution. On the other hand, the models i2.2\ and ( I2.3l l prevent the appearance of 
high frequencies by nonlinear self-interaction. The above mechanism is essentially the one 
that leads to the norm inflation phenomenon in ||4][l4]|22l, except that in those papers, the 
approximation by an ODE is used only on a time interval where the _ff "-norm becomes 
unbounded, but not the iJ'^-norm for any k < s. The above mechanism at the PDE level 
leads to the loss of regularity l[Tl[8l|23]|26l, where indeed k is allowed to be smaller than 
s, as recalled in the introduction. Roughly speaking, the appearance of oscillations is quite 
similar to the above ODE example; in the PDE case, the numerology is different, and the 
proof is more intricate. 

2.2. Choice at the PDE level. We consider now the equations 

(2.4) idtU + P{D)u^eU{\u\^'')u, 
and 

(2.5) idtu + P{D)u = e{U{\u\'^)y u, 

where P{D) is a Fourier multiplier with a real- valued symbol P : H'^ — > R, 

p{D)f = piom- 

The L^-norm of u is formally independent of time: 

d f 

(2.6) — / \u{t,x)\^dx ^0. 

at J^d 

In view of this conservation and of Young inequality 

(2.7) mf)h^ ^WKh^Wfh., 

the option (12.5b seems more interesting than i2.4i . and we have the following result. 

Theorem 2.1. Lef ct e N, e e {±1}, P : R."^ ^ R and x e 5(R'^) even and real-valued. 

• For any uq £ L'^{'R'^), ( 12.51 ) has a unique solution u € C(R; Lp'iRf')) such that 
u\t=o — Uq. Its L^-norm is independent of time: i2.6i holds. 

• If in addition uq G H'i'R'^), s e N, then u e C(R; H'{R'^)). 

• The flow map uq i-^ u is uniformly continuous from the balls in L^(R'') to 
C(R; _L^(R'')). More precisely, for all Uq, Vq G L^(R''), there exists C depending 
on a, \\K\\l<x, ||uq||^2 an^^ |ji;o||i2 such that for all T > Q, 

(2-8) ||u - w||L-([-T,T];L2(Rrf)) \\uo - ^^o||L2(R<i)e'^^, 

where u and v denote the solutions to (12.5b with initial data uq and Vq, respectively. 

• More generally, let s G N. For all uq, Vq G i7'*(R'*), there exists C depending on 

a, \K\\v/s,a^, and \\Vq\}1!> such that for all T > ^, 

(2.9) ||m - 'y||L~([-T,T];ff-^(R<i)) ^ \Wa - fo||H=(Rd)e'^^. 

Remark 2.2. As pointed out in lfT2l . even if the solution is constructed by a fixed point 
argument, the continuity of the flow map is not trivial in general. In the case of Schrodinger 
equations dl.lb . continuity of the flow map in i7''(R'^) is known only in a limited number 
of cases: see |E1 for s = 0, HH for s = 1 and s = 2, and lH!) for < s < 1. 
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Proof. First, recall that S'(t) — e **^(^) is a unitary group on 7? '*(R'*), s £ R. Duhamel's 
formula associated to (12.5b reads 

(2.10) u{t) = S{t)uo-ie S{t - t) (^{K * \u\^y {T)dT. 

The local existence in stems from a standard fixed point argument in 

XiT) - {we C([-T,r];L2(R'^)); \\u\\l^^[_t.t],l^-) ^ 2\\uo\\l^}. 
Denote by <^{u){t) the right hand side of (ITTOl l. In view of ^j}, for t e [-T, T], 



^ \\uo\\l^ 

< II"oIIl2 



{K*\ufyu) (t) 



dr 



LftT * \u(t) 



\\uiT)h2dT 



\\u{r)\ 



'dr. 



By choosing T > sufficiently small, we see that X{T) is stable under the action of $. 
Note that in the case of the model ( I2.4l i. the above estimate would have to be adapted, 
forcing us to work in a space smaller than X{T) (L^ regularity in space would no longer 
be sufficient in general). Contraction is established in the same way: 



\Mu){t) - ^v)it)\U2 ^ I \({K^H^Yu){t)-({K*\v\'Yv){t) 



-T 
T 



dr 



L2 



j \^({K*\u\''Y ^{K*\v\ 
{K*\v\y^ {u~v) 



dr 



T 

-T 



dT. 



L2 



Using the estimate \a' 

mn){t)-mmL^<\\K\\i 

+ \\K\ 



^"1 ^ (lar^^ + W'^^)\a - b\, and (O again, we infer 

T 



\U\ 



\2a-l 
L2 



\v\\l%\\u-vh2dT, 



where all the functions inside the integrals are implicitly evaluated at time r. Choosing 
T > possibly smaller, <i> is a contraction on X{T). Note that this small time T depends 
only on a, \\K\\lc^ and ||uo||£2. Since the L'^-norm of u is preserved (see e.g. ifTTI for 
the rigorous justification), the construction of a local solution can be repeated indefinitely, 
hence global existence and uniqueness at the level. 

Global existence in H'^{'R'^) for s e N then follows easily, thanks to the estimate 



{K*\u\'^yu 



< 



E 



5" (K^M'^Yd^u 



|a| + |^|=s 



L2 



< lli^ll 



6 



R. CARLES 



The continuity of the flow map in is obtained by resuming the estimate written to 
establish the contraction of $: for < > 0, 

\\u{t) - v{t)\\L^ ^ \\uo - v4l^ + \\K\\l^ f {\\uf^. + \\vf^.) \\u - vU^dr 

Jo 

^\\uo-vo\\L^ + \\K\\l^{\\uofj. + \\vo\\T2) [ \\u-v\\L^dT, 

Jo 

where we have used the conservation of the L^-norm. Proceeding similarly for t < 0, 
Gronwall lemma then vields ( 12.81 1 for C depending onlv of a. \\K\\l<x, ||uo||l2 and ||uo||l2- 
Finally, i2.9i is obtained in a similar fashion. □ 

Remark 2.3. The proof of continuity of the flow map is extremely easy. This is in sharp 
constrast with the case of the equation without frequency cut-off. In the case of Schrodinger 
equations (P(^) = continuity is more intricate to establish (see 1271 ). and is true 

only for L^-subcritical nonlinearities, a ^ 2/d, from I.14J . 

We note that even for large a, global well-posedness in is available, in sharp contrast 
with the nonlinear Schrodinger equation dl.lt . Even in the focusing case e = — 1, the high 
frequency cut-off prevents finite time blow-up. In ( 12.91 ), consider vq = v = Q and s = 1 
for instance: by comparison with the case of (11.11 ). we see that the constant C necessarily 
depends on K (or equivalently on x), and is unbounded as x converges to the Dirac mass. 
The frequency cut-off 11 removes the instabilities, and prevents finite time blow-up. 

Remark 2.4 (Hamiltonian structure in the cubic case). If cr = 1, ( 12.4b and ( 12.51 1 coincide. 
We have the equivalence 

X even and real-valued <;=^ K even and real-valued. 

This implies that under the assumption of Theorem l2.1l (12.5b has an Hamiltonian structure, 
and the conserved energy is 

H{u) = J u{x)P{D)u{x)dx + ^ jj K[x-y)\u{y)\^\u{x)\'^dxdy. 

3. Convergence in the smooth case 

Suppose that P{D) converges to A and that 11 converges to Id: does the solution to ( 12. 5b 
converge to the solution of NLS? We show that this is the case under suitable assumptions 
on these convergences, at least in the case where the solution to the limiting equation dl.Ib 
is very smooth. In the sequel, the convergence is indexed by ft, e (0, 1]. 

Proposition 3.1. Let ct G N. We assume that P and 11 verify the following properties: 

• There exist a, (3 ^ such that Ph{0 = + O (h" (O'')- 

• XhiO = X i^O' with X G 5(R'^; [0, 1]) even, real-valued, x — ^ on the unit ball. 
Denote by the solution to ( 12.51 ) with P^ and Xh, such that u^^f^^ = i*|t=o- Suppose that 
the solution fo dl.Ib satisfies u G L°° {[0,T]; H'^'^^), for some s > d/2. Then 

\\U - U ||L~([0,T];ifO " 

Example li.l. The above assumption on Ph is satisfied with a = 1 and /3 = 2 in the 
following cases: 
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• PhiO = -■^arctaii(/i|C|2). 

The second example is borrowed from |fT6l|, where this truncated operator appears naturally 
when discretizing the Laplacian for numerical schemes. 

Remark 3.3. In this result, no assumption is needed on the possible decay of x at infinity. 
Proof. Let = u — u^: it satisfies w[J^q = and 

idtw'^ + Ph{D)w'' = e {n^ {\u\^)Yu - e (H, {\u^\'')Y u'^ 

+ {Ph{D) -^)u + e (l^/p- - {Hh {\u\^)Y) u, 



where we have denoted by II/i the Fourier multiplier of symbol Xh- Denote by R^{u) the 
second line, which corresponds to a source term. In view of the assumption on Ph, there 
exists C independent of /i e (0, 1] such that 

\\PHiD)f ~ M\W ^ V/ e 

We also have, by Plancherel formula, 

II (1 - n,o /f«. = / (1 - xmf io'- \m\'d^ 



j\e\>i/h 

j\i\>i/h 



Therefore, 



P'(")I|l~([0,T];//=) < /^'"'"("•''^lkllL==([0,T];ff.+^). 

Now since s > d/2, H" {YV^) is an algebra, and there exists C independent of h such that 

(n, {\u\^)Yu (n, {Wf)Y ^ ^ cWxWl^ {\\u\\% + \\u^\\%) \\u - u'^\h^, 

where the Young inequahty that we have used is not the same as in Section [2] 

l|i^*/l|L^ ^ 

This is essentially the only way to obtain an estimate independent of G (0, 1]. Indeed, 
Iih{f)=Kh*f, with 

Kn{x)= ^ -^"^ 



The result then stems from a bootstrap argument: so long as 

ll"'*l|L-([0,t];ffO ^ 1 + II"I|l~([0,T];H-), 

Gronwall lemma yields 

Therefore, up to choosing h sufficiently small, this estimate is valid up to t = T. □ 

Such a convergence result can be compared to the one proved in ifTSl to prove the 
convergence of numerical approximations. The approach there is a bit different though, 
inasmuch as the frequency cut-off does not affect the nonlinearity (as in ( 12.5b ). but the 
initial data: consider u'' solution to 
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Then in ifTsl . the discrete analogue to Hhu — Vh is proven to be small. Proposition 13. II 
differs from the results in ifTsl on several aspects: 

• The context in ifTSl is discrete. 

• Only the low frequency part of u, HhU, is shown to be well approximated. 

• The regularity assumption on u may be much weaker 

As mentioned above, the second point is due to the choice of the model. However, as 
discussed in Section 12.11 controlling the high frequencies of the initial data must not be 
expected to ensure a control of high frequencies of the solution w'' for positive time. 

The third point is due to the use of Strichartz estimates in ifTSll . In the next section, 
we show that in the presence of dispersion (with Ph{£,) = Proposition 13 . II can be 

adapted to rougher data. 

4. Convergence using dispersive estimates 
We first recall the standard definition. 
Definition 4.1. A pair (p, q) ^ (2, oo) is admissible if p ^ 2, q ^ 2, and 



We shall consider (I2.5l l when P{D) is exactly the Laplacian, and not an approximation 
as in Proposition 13. II The reason is that when P is bounded, then no Strichartz estimate 
is available, as we now recall. Let S{-) be bounded on PP^ for all s ^ 0. By Sobolev 
embedding, for all {p, q) (not necessarily admissible) with 2 ^ q < (x, there exists C > 
such that for all uq E i7''/^~'*/'?(R'^), and all finite time interval /, 

ll'S'(-)"o||LP(/:L<!(Rd)) < C||S'(-)wo||LP(/;//<*/2-d/9(R'i)) 

^ C'll'"0||LP(/;ffd/2-d/,(Rd)-) = C|/|^/^||Uo||//d/2-d/5(Rd). 

If the Fourier multiplier P is bounded, the above estimate cannot be improved, in sharp 
contrast with the result provided by Strichartz estimates. 

Proposition 4.2 (From |9 |). Let d ^ 1, and P e L°°{R'^; R). Denote S{t) = e-'^^'^^l 
Suppose that there exist an admissible pair (p, q), an index fc g R, a time interval I 3 0, 
\I\ > 0, and a constant C > such that 

WSi-hohni-.LH-R''}) ^ C'||Mo||H'=(R<i), Vuo e H'^i'R'^). 
Then necessarily, k ^ 2/p — d/2 — d/q. 

We now state the main result of this section. 

Tlieorem 4.3. Let a eN and T > 0. We assume that Xh{0 = X {^0' ^^^^ X S S{Hf-) 
even, real-valued, x = 1 on B{0, 1). Let u solve (II. lb , and consider the solution m'' to 

1. Suppose that a ^ I and d 2. Ifu £ L°°([0, T]; L^) n L^^'^{[0,T]; L'^), then 

w''||l-([0.T];L2)^— ^0. 

2. Suppose that (7 = 1 and d ~ 3. 

• Ifu,\7u £ L°^{[0,T]; L^) n L^/'^{[0,T]; L% then 

\\u-u^\\L^^[o,T];m)^^0. 



NLS AND FREQUENCY SATURATION 



9 



• IfueL^{[0,T];H'), with s > 3/2, then 

I|w-'«''IIl~([o,t];L2) < /i" and \\u - u^\\l^(^[o,t]-h') < h'^^ ■ 
3. Suppose that cr ^ 1 and d ^ 2. If u £ L°°([0, T]; if"), with s ^ 1 and s > d/2, then 

||m - '"'NIl~([o,t];L2) < ft-^ and ||u - u''||i=c([o,T];ifi) ^0. 

If in addition s > 1, then 

||u-u'*||ioo([o,r];Hi) < h''-^. 

Remark 4.4. Suppose uq sufficiently smooth. If e = +1 (defocusing case), the bounds for 
u are known in several cases, with T > arbitrarily large. On the contrary, if e — — 1 
(focusing case), T may have to be finite, bounded by a blow-up time. See e.g. ifTTlfrTl . 
Typically, if cr = d = 1, then the assumption of the first point is fulfilled for all T > as 
soon as uq G L^(R), for e G {±1}, from |27|, and if cr 1, ^ 2, the assumption of the 
third point is fulfilled for all T > as soon as uq G i/'' (R'*), for e = +1, from fTTI . 

Proof For fixed h > 0, Theorem lO shows that u'^ G C{'R;H''), with /c = 0, 1 or 
s according to the cases considered in the assumptions of the theorem. Of course, the 
bounds provided by Theorem l2. ll blow up as /i — > if fc > 0. 

As in the proof of Proposition 13. II let = u — u'\ The equation satisfied by is 
simpler than in the proof of Proposition l3.1l since Ph{D) = A: 

^^tw'^ + AuA ^ e (n, (1^.^)% - e (H, (1^'^^)%'^ + 6 [\u\"^ - (H, {\u\')y) u. 

Resume the notations R'^{u) = e (lup"" - {Uh u and Ilh{f) = Kh * /, with 

Kh{x) = {2iT)-'^/^h-'^x{~x/h). From young inequality, we have, for all g G [1, oo], 

(4.1) l|n„(/)|U, ^ ||A-,|Ui||/|U, ^ \\x\\LA\f\\L., 

an estimate which is uniform in /i > 0. Introduce the Lebesgue exponents 

4(7 + 4 2CT(2cr + 2) 



9 = 2cr + 2 ; p 



da ' 2-{d~2)cr' 



The pair (p, q) is admissible, and 

1 2f7 1 1 2cr 1 

(4.2) ^_ + _ ; = +_. 

q q q p' t) p 

For t > 0, denote L^L'' = U{[0, t]; L'=(R'')). From Strichartz estimates (see e.g. IfTlT ). 



< 



(ll«lli«L, + ll"'lli«L,)lk'llLrL, + ||i?'N| 



where we have used Holder inequality and (14. Il l, and where (pi, gi) is an admissible pair 
whose value will be given later. 

\f a — 1 and d ^ 2, then 6 ^ p, and we infer 



In the first case of the theorem, we assume u G LP{[0,T]; L'^), since p — 8/d and q ~ 4 
for (7 = 1. We use again a bootsti-ap argument: so long as IIu'^Hlpl? ^ 2||u||iPi<j, we 
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divide the interval [0, T] into finitely many small intervals so the first term of the right hand 
side is absorbed by the left hand side (recall that p is finite), and we have 



The bootstrap argument is validated provided that II II p/ , — as /i — > 0. 

If we have only a < -^^j then by Sobolev embedding. 

In the same way as above, 

l|v«^'^L.i,nLrL^ ^ II V ((n, (|«|^))% - (n, {\u'^f)Y u'^ 

II ^ ^IlLfiL'i 

The first term of the right hand side is controlled by 

{nu{\u\^)YVu-{ll,,{\u^\^)Yyu'^ 



L^J Li' 



(4.3) 



Introducing the factor (ll/i (jup))'^ Vu'*, the first term is estimated by 



< 



< 



re/2 /2 



Li' Li' 



{n,{\uf)y-{n,{\u'^f)y) vu^ 



Li' Li' 



L1L1 



{\\<rLi + W^XTl") III"!' - l^'HLrL,/^ 



Li 



L»Li 



\Vw''\ 



LiLi 



2(7-1 

Li 



hu2a-l 
^ "LfLi 



KL^Li W^^ WlILi 



\^w''\\l?Li 



+ t 



2a 16 



\2u-l 



/i||2ct-1 



LiLi 



Proceeding similarly for the other term in (14.3b . splitting [0, T] into finitely many time 
intervals where the terms containing on the right hand side can be absorbed by the left 
hand side, and using a bootstrap argument, we end up with 

\\w''\\Liw^.inLTm < . 

Therefore, it suffices to show that for some admissible pair (pi, qi), the source term con- 
verges to in iP'i ([0, T]; L«i) (if ct = 1 and d < 2) or in L^'i ([0, T];W'^'i'^) (in the other 
cases), so the bootstrap argument is completed. In addition, the rate of converge of the 
source term, if any, yields a rate of convergence for . The theorem then stems from the 
following lemma, in which (p, q) is given by ( I4.2l l. 

Lemma 4.5. Let T > 0. The source term R^{u) can be controlled as follows. 

1. Suppose that a = I and d 2. Ifu £ L°°([0, T]; L^) n L*^/'*([0, T]; L"*), then 

II-K''(w)IIlp'([o,t1;L<!') ^^0- 

2. Suppose that a ~ 1 and d — 3. 
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• Ifu,Vu e L°^{[Q,T];L^) rM^'<^{[Q,T]]L^), then 

II^'*(")IIlp'([0,T];W1'9') ^q'^- 

• Ifu e i°°([0,r];iJ''), with s > 3/2, then 

\\R'\u)\\m[o.T];L-) < and ||i?"(u)|Ui([o,T];ffi) < 
3. Suppose that a ^ 1 and d ^ 2. Ifu £ L°°([0, T]; if"), with s ^ 1 ant/ s > d/2, then 

If in addition s > 1, then 

Proof of Lemma \43\ For the first case, we use Holder inequaUty, in view of i4.2V . 

= IKi-n/.)(l"l')«IL.'i,' ^ 11(1 -n/.)(kP) 11^/2^,/. Iklli.L,- 

We note that for cr = 1, g = 4, so by Plancherel Theorem, 

II (1 - n,) (i^^niii, = / (1 - xmf \Ti\u\^)m'd^ < / \H\nm)?dt 

By assumption, M G LP{[0,T]; L'^) C L^([0, T]; L^)^ thus |u|2 e L^/^dO, T]; L^), and by 
Plancherel Theorem, J^(|Mp) € I/^/^([0, T]; L^). The first point of the lemma then stems 
from the Dominated Convergence Theorem. 

For the first case of the second point, we note that now 9 > p, so the above argument 
must be adapted, and we have to estimate the gradient of R^{u) in the same space as above. 
Since L°°([0, T]-H^[Ii^)) C L%[0, T]]L^{Ii^)), the Dominated Convergence Theorem 
yields 

We now estimate WR'^{u). Write 

+ ||(i-n/,)v(|up)||^(i/,+i/,,-i^^ hllL«L2 
+ 11(1- n,o V {\uf) Ml^l^ 

By the same argument as above, 

||(i-n;,)(|^.p)||^^^, IIVuLp^.^^o. 

We note that u bounded in L°°([0, T];H\R^)) C L^([0, T]; L'^(R^)), and Vu bounded 

in L^L'^, so V|mP is bounded in l'^^^^^^^'' L^. Invoking Plancherel Theorem and the 
Dominated Convergence Theorem like above, we infer 

||(i-n,)v(|zip)||^(,,«^,/,,-i^^ l|y||^^^,_,o. 

This completes the proof for the first case of the second point. 
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For the remaining cases, we use that H^H'^) is embedded into L°°(R'*): for fixed t. 



P'HWIU^ < {\Ht)\\ 



2cr-2 



|n.(Kt)nill;.^)||(i-n,)(Ki)niU.||^(t)lh 



< ll"Wlir=o"'ll(i-n,)(Kt)niU. <||u(i)f/r^||(i-n,)(K<)p^ 

Like in the proof of Proposition l3.1l we use the estimate 
(4.4) ||(l-n,)/|U2 s^/i^ll/llff., 

and since i?'*(R'*) is an algebra. 



To conclude the proof, we estimate VR'''{u) in L^(R''). We compute 

\7R\u) = a|up--2 _ (V u 

+ a[\ur-'-{n,{\u\')y-')n,{v{\u\'))u 

+ [\u\''^-{U,{\u\')y) Vu, 
where the second line is zero is ct = 1. We estimate successively, thaks to ( 14. It . 

II ((1 - n,o (V (lup))) u||^, < \\u\\l^-' 11(1 - n,o {\u\') \\^, , 

I (l^p- - (n, vz^ll^^ < \\u\\l^-' 11(1 - n,o {\u\^) 11^^ iiv^iu., 

and, if (T ^ 2, 



(H^-^-(n.(|.p))^-^)n. (v(|z.n)^ 



L2 



<||zi|li'L-^||(i-n,)(|«n 



|(i-n„)(|«| 

Since we have _ff*(R'') ^ L°°(R''), we end up with 

\\yR\u)h.<\\u\\%-'\\{l~U,){\u\^)\ 
If s > 1, (14.41 ) yields, since in addition s > d/2. 



11(1 -n,) (1^.1^ 



\u\\h 



If s = 1 (a case which may occur only if d = 1, since s > d/2), we write 



|v(i-n,)(|«p)||^, ^ 



\:F{v{\u\^)){0\'dt 



Nowsince V = 2ReMVMandM e H'^i'R) ^ L°°(R), Vu e L2(R), we conclude 

thanks to the Dominated Convergence Theorem. □ 



This completes the proof of Theorem l4.3l by choosing {pi,qi) — {p, q) or (oo, 2). □ 
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Appendix A. Physical saturation of the nonlinearity 

Instead of cutting off the high frequencies, one may be tempted by saturating the non- 
linear potential, by replacing |up not by n(|up) but by /(|itp) where / is smooth, equal 
to the identity near the origin, and constant at infinity. Note also that a saturated nonlin- 
earity may be in better agreement with physical models (recall however that (II. lb appears 
in rather different physical contexts, such as quantum mechanics, optics, and even fluid 
mechanics), since typically the power-like nonlinearity in ( II. lb may stem from a Taylor 
expansion; see e.g. 12111251 . More precisely, let / G C°°(R; R) such that 

sif s; s 1, 

2 if s ^ 2. 



(A.l) f{s) = 



h{\un = lf{h\u\') 



The analogue of the Fourier multiplier H/i is defined as 
and we replace (12.51) with 

(A.2) idtu'' + Ph{D)u'' ^e{h{\u''fyu\ 

so the formal conservation of the L^-norm still holds. We could also consider 

(A.3) A(|„H = _J^!!_^. 

In both cases, the main aspect to notice is that fh is bounded and z H> fh{\z\'^)°'z is 
globally Lipschitzean. We infer the analogue of Theorem l2.1l at least in the case. 

Proposition A.l. Let treN, eG{±l}, P: R'^ R and f given either by (lA.lb or by 
(lA3t . 

• For any Uq G L'^iJEif'), (I A. 2b has a unique solution G C(R; L^(R'')) such that 
it|*j^P = uq. Its L^-norm is independent of time. 

• The flow map uq i—^ is uniformly continuous from the balls in L^(R'^) to 
C(R; L^(R'')). More precisely, for all uq, Wq G i^(R''), there exists C depending 
on a, h, 1 1 Wo 1 1 l2 and ||wo 11^2 such that for all T > 0, 

- t>''||ioo([_7^^7^].i2(R,d)) < ||uo - 'yo||L2(Rd)e'^'^, 

where and denote the solutions to (IA.2b with data uq and vq, respectively. 

Introduce 

Fhis) = I IhivYdy. 
Jo 

We check that the following conservation of energy holds: 

u'\t,x)Ph{D)u''{t,x)dx + e J ^Fh{\u{t,x)\^)dx^ =0. 

Proving the analogue of Proposition l3.1l is easy in the case (lA.lb . since the last source term 
for the error w'' is now 

R'\u)={\u\''^-f,{\u\'y)u, 

and under the assumptions of Proposition |3.1| u G L°° ( [0, T] x R''), so there exists ho > 
such that for < h ^ ho, 

\u{t, = {\u{t, x)\^y , V(<, x) G [0, T] X R'^. 
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Therefore, this source term simply vanishes for h sufficiently small. In the case (IA.3b . we 
can use the relation 



(A.4) \Ru{u)\^\{\u\^'^ -U{\ufy) 

and Schauder lemma to get a source term which is 0{h) in i7'*(R'*), for s > d/2. 

Proposition A.2. Lef cr G N. We assume that P is such that there exist a, /3 ^ with 
Ph{C) = + O {h'^ iO^y Denote by the solution to (IA.2b w;f/z Ph and fh, such 

that uj'j^g — U|j^o- Suppose that the solution to jl.ll l satisfies u £ L°°([0, T]; H'^^^), for 
some s > (i/2. 

• /« f/ze ca.se ( lA.ll l, ||w — 'u^\\l°°{[o.t]-H'') ^ 

• In the case (A3, \\u - u''\\l^(^[o,t]-H') < 

In the case jA.ll ). proving an analogue to Theorem |4.3| seems to be more delicate though, 
and we choose not to investigate this aspect here. On the other hand, in the case (IA.3I ), 
using the estimate ( IA.4b . Strichartz estimates and Holder inequalities with the "standard" 
Lebesgue exponents (in the same fashion as in the proof of Theoi'em l4.3l see e.g. lillJ ). we 
have, with steps similar to those presented in the proof of Theorem l4.3l 

Tlieorem A.3. Let ct e N and T > 0. Let u solve ( 11.11 ), and consider solution to 

1. If a ^ 2/d, andu£ L°°{[0,T]; L^) D 1^"^"+'^^/'^'' {[0,T]; L^''+'^), then 

I|w-w''IIl-([o.t];L2)^-^0. 

2. Suppose that a = 1 and d = 3. 

• //u, Vu e L°^{[Q,T]\L^) C\ L^/'^{[Q,T];L% then 

• Ifu e i°°([0, T]; iJ"), with s > 3/2, then 

3. Suppose that a ^ 1 and d ^ 2. Ifu £ L°°{[0, T]; H'^), with s ^ 1 and s > d/2, then 
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